o 

(N 
>^ 



o\ 






> 

m 

rn 

in 

o 
o 



X 



Single-shot qubit readout in circuit Quantum Electrodynamics 

Frangois ^Mallet, Florian R. ^Ong, Agustin ^Palacios-Laloy, Frangois 

^Nguyen, Patrice ^Bertet, Denis ^Vion and Daniel ^Esteve 

^Quantronics group, Service de Physique de I'Etat Condense (CNRS URA 2464), 

DSM/IRAMIS/SPEC, CEA-Saclay 
91191 Gif-sur-Yvette Cedex, Francq^ 



"Electronic address: ldenis.vion@cea.f'rl 



The future development of quantum information using superconducting cir- 
cuits requires Josephson qubits [l] with long coherence times combined to a 
high-fidelity readout. Major progress in the control of coherence has recently 
been achieved using circuit quantum electrodynamics (cQED) architectures (2,|3|, 
where the qubit is embedded in a coplanar waveguide resonator (CPWR) which 
both provides a well controlled electromagnetic environment and serves as qubit 
readout. In particular a new qubit design, the transmon, yields reproducibly 
long coherence times jj, |5|. However, a high-fidelity single-shot readout of the 
transmon, highly desirable for running simple quantum algorithms or measur- 
ing quantum correlations in multi-qubit experiments, is still lacking. In this 
work, we demonstrate a new transmon circuit where the CPWR is turned into 
a sample-and-hold detector, namely a Josephson Bifurcation Amplifer (JBA) 
(g, |7J , which allows both fast measurement and single-shot discrimination of the 
qubit states. We report Rabi oscillations with a high visibility of 94% together 
with dephasing and relaxation times longer than 0.5 jis. By performing two 
subsequent measurements, we also demonstrate that this new readout does not 
induce extra qubit relaxation. 

A common strategy to readout a qubit consists in coupling it dispersively to a resonator, 
so that the qubit states |0) and |1) shift differently the resonance frequency. This frequency 
change can be detected by measuring the phase of a microwave pulse reflected on (or trans- 



mitted through) the resonator. Such a met 
pair box capacitively coupled to a CPWR 



lod, successfully demonstrated with a Cooper 



Ql 



faces two related difficulties which have 



prevented so far from measuring the qubit state in a single readout pulse (so-called single- 
shot regime): the readout has to be completed in a time much shorter than the time Ti in 
which the qubit relaxes from |1) to |0), and with a power low enough to avoid spurious qubit 
transitions J8|. 

This issue can be solved by using a sample-and-hold detector consisting of a bistable 
hysteretic system whose two states are brought in correspondence with the two qubit states. 
Such a strategy has been implemented in various qubit readouts 



% 



a 



n our experi- 



ment the bistable system is a Josephson Bifurcation Amplifier (JBA) p, |7| obtained by 
inserting a Josephson junction in the middle of the CPWR (see Fig. 1). When driven by 
a microwave signal of properly chosen frequency and power, this non-linear resonator can 



bifurcate between two dynamical states B and B with different intra-cavity field amplitudes 
and reflected phases. In order to exploit the hysteretic character of this process, we perform 
the readout in two steps (see Fig. [1^): the qubit state |0) or |1) is first mapped onto B or 
i? in a time much shorter than Ti, the selected resonator state is then hold by reducing the 
measuring power during a time tn long enough to determine this stateQ with certainty. 

JBAs were used previously to readout quantroniums ll|-|l3| and fiux-qubits, obtaining 
for the latter fidelities up to 87% jlj] with Quantum- Non-Demolition character 15|. Here 



we couple capacitively a transmon to a JBA, combining all the advantages of the cQED 
architecture (long coherence times, scalability) with the single-shot capability of a sample- 
and-hold detector. A crucial characteristic of this new design is its very low back-action 
during readout. Indeed the qubit frequency depends only on the slowly- varying photon 



number inside the resonator 16J, yielding less relaxation than in previous experiments where 
the qubit was coupled to a rapidly varying variable of the JBA (the intra-resonator current). 
Furthermore we designed the resonator to make it bifurcate at a low photon number, thus 
avoiding unwanted qubit state transitions during readout. 

The complete setup is shown in Fig. [1] the transmon |4,|5j of frequency /oi tunable with a 
magnetic flux (p is coupled with a coupling constant g = 44±3 MHz to the non-linear CPWR 
of fundamental frequency /c = 6.4535 GHz, quality factor Qq = 685 ±15 and Josephson 
junction critical current Iq = 0.72 ± 0.04 pA. In this work the qubit is operated at positive 
detunings A = /c — /oi larger than g. In this dispersive regime the resonator frequency /ci 
depends on the qubit state \i), and the difference 2% = fco — fci (so-called cavity pull) is 
a decreasing function of A. Readout pulses (Fig. [1^) of frequency / and maximum power 
Ps are sent to the circuit; after reflection on the resonator their two quadratures I and 
Q are measured by homodyne detection. They belong to two clearly resolved families of 
trajectories (Fig. [Tb) corresponding to both oscillator states B and B. The escape from 
S to S is a stochastic process activated by thermal and quantum noise in the resonator 

nn 

|17l . Il8| . and occurs during the sampling time ts with a probability ps that increases with 
Ps- The position of the so-called "S-curve" pb(-Ps) depends on the detuning /ci — / (6| and 
thus on the qubit state. When the two S-curves Sf and S^ corresponding to |0) and |1) are 
sufficiently separated, one can choose a value of Ps at which these states are well mapped 
onto B and B (Fig. [It). 

We now present our best visibity, obtained at A = 0.38 GHz in this work and confirmed 



on another sample. We measure 5*^ and S^ (Fig. [2]) after preparing the transmon in state |0) 
or |1) using a resonant microwave pulse. The contrast, defined as the maximum difference 
between both curves, reaches 86%. To interpret the power separation between the S-curves, 
we search the readout frequency /+ Afi that makes 5'f'^^f^ coincide with S^ at low bifurcation 
probability. This indirect determination of the cavity pull gives Afi = 4.1MHz, in good 
agreement with the value 2% = 4.35 MHz calculated from the experimental parameters. At 
high pb however the two S-curves do not coincide, which reveals that the limiting factor of 
our readout fidelity is relaxation of the qubit before the time needed for the resonator to 
reach its final state. To reduce this effect and improve the readout contrast, we transfer 
state |1) into the next excited state |2) with a resonant vr pulse just before the readout 
pulse, yielding the S-curve Sf and a 92% contrast. This technique, already used with other 



Josephson qubits lOJ, is analogous to electron shelving in atomic physics and relies here on 
the very low decay rate from |2) to |0) in the transmon. Figure [2}3 shows Rabi oscillations 
between |0) and |1) obtained with such a composite readout pulse. The visibility, defined 
as the fitted amplitude of the oscillations, is 94%, and the Rabi decay time is 0.5 ]as. Of 
the remaining 6% loss of visibility we estimate that about 4% is due to relaxation before 
bifurcation and 2% to residual out-of-equilibrium population of |1) and to control pulse 
imperfections. Such a visibility higher than 90% is in agreement with the width of the S- 
curves estimated from numerical simulations, with their theoretical displacement, and with 
the measured qubit relaxation time. 

The visibility being limited by relaxation, it is important to determine whether the read- 
out process itself increases the qubit relaxation rate. For that purpose we compare (at 
A = 0.25 GHz) Rabi oscillations obtained with two different protocols: the control pulse 
is either followed by two successive readout pulses yielding curves Ri and R2, or by only 
the second readout pulse yielding curve R3 (see Fig. [3^). R2 and R3 exhibit almost the 
same loss of visibility compared to Ri, indicating that relaxation in the presence of the first 
readout pulse is the same as (and even slightly lower than) in its absence. 

To further investigate this remarkable effect, we measure Ti in presence of a microwave 
field at the same frequency / as during readout, and for different input powers P (see Fig. 
|3h>). We first roughly estimate the intra-cayity mean photon number n{P) by measuring 



the AC-Stark shifted qubit frequency /oi(-P) 16| (the correpondence fm{n) is obtained by a 



numerical diagonalization of the Hamiltonian of the transmon coupled to a field mode with 



n photons). Bifurcation is clearly revealed by a sudden jump of n from about 5-10 to 50-100 
photons. Meanwhile Ti does not show any decrease up to about 5 dB above bifurcation. It 
even slightly increases because the qubit frequency is pushed away from the cavity, slowing 
down spontaneous emission as explained in the next paragraph. This is in strong contrast 
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19|. These results prove that our 



with all previous experiments using a JBA readout 
design achieves very low back-action on the qubit. A similar behavior was observed for most 
qubit frequencies, except at certain values of P and /oi where dips in Ti (P) were occasionally 
observed above bifurcation. 

We now discuss the dependence of the readout contrast and qubit coherence on the de- 
tuning A. Besides acting as a qubit state detector, the resonator serves also as a filter 
protecting the qubit against spontaneous emission into the 50 Vt impedance of the exter- 
nal circuit 20|, |21|. The smaller A, the stronger the coupling between the qubit and the 
resonator, implying a larger separation between the S^ and S} curves but also a faster re- 
laxation. We thus expect the contrast to be limited by relaxation at small A, by the poor 
separation between the S-curves at large A, and to exhibit a maximum in between. Figure 
m presents a summary of our measurements of contrast and coherence times. At small A, 
Ti is in quantitative agreement with calculations of the spontaneous emission trough the 
resonator. However it presents a saturation, similarly as observed in previous experiments 



20|, but at a smaller value around 0.7 ps. The effective cavity pull Afi determined from the 
S-curves shifts (cf. Fig. [2]) is in quantitative agreement with the value of 2% calculated from 
the sample parameters. The contrast varies with A as anticipated and shows a maximum 
of 92% at A = 0.38 GHz, where Ti = 0.5 ps. Larger Ti can be obtained at the expense of 
a lower contrast and reciprocally. Another important figure of merit is the pure dephasing 
time T^ 23| which controls the lifetime of a superposition of qubit states. T^ is extracted 
from Ramsey fringes experiments (see Methods), and shows a smooth dependence on the 
qubit frequency, in qualitative agreement with the dephasing time deduced from a 1// flux 



noise of spectral density set to 2O]a0o/vHz at 1 Hz, a value similar to those reported else- 



where [2^. To summarize our circuit performances, we obtained a 400 MHz frequency range 
(pink area on Fig. H]) where the readout contrast is higher than 85%, Ti is between 0.7 ps 
and 0.3 ps, and T^ between 0.7 ps and 1.5 ps. Further optimization of the JBA parameters 
Jc and Qq could increase this high-visibility readout frequency window. 

In conclusion we have demonstrated the high-fidelity single-shot readout of a transmon 



qubit in a circuit-QED architecture using a bifurcation amplifier. This readout does not 
induce extra qubit relaxation and preserves the good coherence properties of the transmon. 
The high fidelity achieved should allow a test of Bell's inequalities using two coupled trans- 
nions, each one with its own JBA single-shot readout. Moreover, our method could be 
used in a scalable quantum processor architecture, in which several transmon- JBAs with 
staggered frequencies are read by frequency multiplexing. 

Methods 

A. Sample fabrication 

The sample was fabricated using standard lithography techniques. In a first step, a 
120nm-thick niobium film is sputtered on an oxidized high-resistivity silicon chip. It is 
patterned by optical lithography and reactive ion etching of the niobium to form the coplanar 
waveguide resonator. The transmon and the Josephson junction of the CJBA are then 
patterned by e-beam lithography and double-angle evaporation of two aluminum thin-films, 
the first one being oxidized to form the junction tunnel barrier. The chip is glued on 
and wire-bonded to a microwave printed-circuit board enclosed in a copper box, which is 
thermally anchored to the mixing chamber of a dilution refrigerator at typically 20 mK. 

B. Electrical lines and signals 

Qubit control and readout microwave pulses are generated by mixing the output of a 
microwave source with "DC" pulses generated by arbitrary waveform generators, using DC 
coupled mixers. They are then sent to the input microwave line that includes bandpass 
filters and attenuators at various temperatures. The powers given in dB in this letter are 
arbitrarily refered to 1 mW (on 50 Q) at the input of the dilution refrigerator; the total 
attenuation down to the sample is about — 77dB. The pulses are routed to the resonator 
through a circulator to separate the input and output waves. 

The readout output line includes a bandpass filter (4 — 8 GHz), 2 isolators, and a cryogenic 
amplifier (CITCRYO 1-12 from California Institute for Technology) with 38 dB gain and 
noise temperature Tn = 3 K. The output signal is further amplified at room-temperature 
with a total gain of 56 dB, and finally mixed down using an I/Q mixer with a synchronized 
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local oscillator at the same frequency. The I and Q quadratures are further amplified by 

20 dB, and sampled by a fast digitizer. The data are then transferred to a computer and 
processed. The single-shot traces of Fig. [T]d. were obtained with an additional 10 MHz 
low-pass filter. 

C. Sample characterization 

The characteristic energies of the system, namely the transmon Josephson energy i?j = 

21 GHz and charging energy Ef. = 1.2 GHz (for a Cooper-pair), as well as the qubit-resonator 
coupling constant g , have been determined by spectroscopic measurements. The bare res- 
onator frequency fc was determined at a magnetic field such that the qubit was far detuned 
from the resonator. 

D. Qubit state preparation 

We prepare the qubit in its ground state with a high fidelity at the beginning of each 
experimental sequence by letting it relax during about 20iis. We estimate at about 1% the 
equilibrium population in state |1) due to residual noise coming from measurement lines. 

To prepare the qubit in its excited state |1) or |2), one or two successive resonant square- 
shaped pulses of length t-j^ ~ 20 ns are applied prior to the readout pulse. The dotted blue 
S-curve of Fig. [T]was recorded with a single resonant vr pulse at /12 (see text): it reveals 
that this pulse induces a spurious population of the |1) state of order 1%. We checked that 
this effect is corrected by using gaussian-shaped pulses J9| (data not shown). 

E. Readout Pulses 

We give here more information on the timing of the readout pulses used is this work. 
In Fig. O readout is performed at /c — / = 17 MHz, and we used tn = 15 ns, ts = 250 ns 
and tn = 700 ns. We stress that although tg is of the same order of magnitude as Ti, the 
observed relaxation-induced loss of contrast is rather low, which may seem surprising. This 
is due to an interesting property of our readout : when the qubit is in state |1), the JBA 
bifurcates with a high probability, implying that all bifurcation events occur at the very 



beginning of the readout pulse (instead of being distributed exponentially during ts). We 
nevertheless keep ts = 250 ns because the bifurcation process itself needs such a duration to 
develop properly. The effective measurement time t^ is thus shorter than ts- We verified 
that weighted sums of Sf and iS'°_,_^f. fit properly the S^ curves (i=l,2) of Fig. [21 allowing us 
to quantify the population of each level at readout. Using the experimentally determined 
relaxation times Tf~^^ ~ 0.3 pis and T^^^ ~ 0.45 ps, we thus estimate tM ~ 40ns. 

In Fig. 121 readout is performed at /c — / = 25 MHz, to reduce the total measurement 
duration. Indeed, as a larger readout detuning implies a higher driving power and thus 
a higher reflected power, the signal to noise ratio is increased which allows to shorten tn 
to 50 ns. We also used for these data t^ = 10 ns and ts = 40 ns to shorten the overall 
measurement time, which also decreases the maximal contrast to approx 83%. Finally, a 
delay time of 120 ns between the two readout pulses has been optimized experimentally 
to empty the resonator of all photons due to the first measurement, and thus avoid any 
spurious correlations between the two outcomes of the sequence. 

F. Coherence times measurement 



The qubit coherence times are measured using standard experimental sequences [22|. For 
the relaxation time Ti, we apply a vr pulse and measure the qubit state after a variable delay, 
yielding an exponentially decaying curve whose time constant is Ti. The coherence time T2 
is obtained by a Ramsey experiment: two 7r/2 pulses are applied at a frequency slightly 
off-resonance with the qubit and with a variable delay; this yields an exponentially damped 



oscillation whose time constant is T2. We then extract the pure dephasing contribution T^ to 
the quantum coherence (a 



the quantum coherence (as well as the associated maximum uncertainty) using the relation 
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Figure 1: Principle of a single-shot readout for a transmon qubit. A transmon (magenta) 
is capacitively coupled to a coplanar resonator (green grayed strips) made anliarmonic by inserting 
a Josephson junction (green cross) at its center. This qubit is coherently driven by a source Vq and 
measured by operating the resonator as a cavity JBA: a microwave pulse with properly adjusted 
frequency / and time dependent amplitude (rise, sampling, and holding times tR, isj and t^, 
respectively - see inset a and Methods) is applied by a second source Vr; this pulse is reflected by 
the system and routed to a cryogenic amplifier and to a homodyne detection circuit yielding the 
two quadratures / and Q. During the "sampling" time ts the electromagnetic field in the resonator 
has a probability pB to bifurcate from a low amplitude state ^ to a high amplitude one B, both 
states corresponding to different amplitudes of / and Q. The "holding" time in is then used to 
average I{t) and to determine with certainty if the resonator has bifurcated or not. (b) Oscillogram 
showing filtered /(i) traces of both types (obtained here with t/j = 30 ns and tg = tu = 250 ns). 
(c) The probability pB depends on / and on the sampling power Ps. The two qubit states |0) 
and |1) shift the resonator frequency, resulting in two displaced S-curves 5'^ and S^ . When their 
separation is large enough, Ps can be chosen (vertical dotted line) so that B and B map |0) and 
|1) with a high fidelity. 
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Figure 2: Best single-shot visibity obtained at A = 0.38 GHz and f^ - f = 17 MHz. (a) 

S-curves Pb{Ps) obtained with the qubit prepared in state |0), |1), or |2) (solid lines S^, Sj and Sf, 
respectively) with the proper resonant vr pulses (top diagram). The maximum differences between 
Sf and 5f (red vertical line) and between the Sf and 5^ (green vertical line) define two readout 
contrasts of 86% and 92%. The readout fidelity is thus increased by using a composite readout 
where the measurement pulse is preceded by a vr pulse at frequency /12 that transfers |1) to |2). 
The dotted blue curve obtained after a single vr pulse at frequency /12, starting from |0), shows 
that this technique has almost no effect on |0). Also plotted are the curves obtained for |0) when 
shifting the readout frequency / by Af^ = 4.1 it 0.1 MHz (red dashed line) and Af2 = 5.1 it 0.1 
MHz (green dashed line) in order to match at low pB the curves obtained for |1) and |2). The 
difference between the corresponding solid and dashed curves is a loss of visibility mostly due to 
qubit relaxation before bifurcation, (b) Rabi oscillations at 29 MHz measured with the composite 
readout, as sketched on top. Dots are experimental values oi p^i^t) whereas the solid line is a fit 
by an exponentially damped sine curve with a 0.5 ps decay time and an amplitude of 94% (best 
visibility). The total errors in the preparation and readout of the |0) and |1) states are 2% and 
6.5% respectively. 
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Figure 3: Effect of the readout process onto the qubit at A = 0.25 GHz and 

fc — f = 25 MHz. (a) Rabi oscillations p^i^t) obtained at Ps = — 30.5dB with the protocols 
sketched on top, i.e. with two successive readout pulses placed immediately after the control Rabi 
pulse (red and blue dots), or with the second pulse only (green dots). The loss of Rabi visibility 
between the red curve (83%) and the blue (44%) and green (37%) ones is due to qubit relaxation 
during the first readout or the delay, (b) Top panel: Spectroscopic determination of the qubit 
frequency /oi when it is AC-Stark shifted by an auxiliary microwave with frequency / and power 
P (protocol on top). The shift provides an in-situ estimate of the average photon number n in 
the resonator (right scale) with a precision of ±30%. The bifurcation is seen as a sudden jump. 
Bottom panel: qubit relaxation time Ti (measurement protocol not shown) in presence of the same 
auxiliary field. Ti does not show any strong defc8ease even at power well above bifurcation. 




foi (GHz) 



Figure 4: Trade-off between qubit coherence and readout fidelity, (a) Experimental re- 
laxation time Ti (red dots) and dephasing time T,^ (violet dots) of the qubit as a function of /oi 



231 (A w -0.75 GHz, 



(or equivalently A/g). Note that T^ ~ 2.5 it 0.5 jis at the flux optimal point 
data not shown). Error bars on T^ are absolute minima and maxima resulting from the maximum 
experimental uncertainties on the coherence times Ti and T2 (see methods). The solid red line is 
the value of Ti obtained by adding to the expected spontaneous emission through the resonator 
(dashed red line) a relaxation channel of unknown origin with Ti = 0.7 jis (horizontal dotted line). 
The blue line is the pure dephasing time T^ corresponding to a 1// flux noise with an amplitude 
set to 20]J(/)o/vHz at 1 Hz. (b - left scale) Readout contrast with (green dots) and without (green 
circles) transfer from state |1) to |2) (see Fig. 2). (b - right scale) Effective cavity pull Afi (blue 
squares) determined as shown in Fig. [2l For the sake of comparison, the predicted cavity pull 2% 
in the dispersive approximation is also shown as a cyan region, taking into account the maximal 
experimental uncertainty on g. The pink area denotes the region where the readout contrast is 
higher than 85%. 
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